P=[2,2];
[P,iter,err]=Newdim('F','JF',P,1.0e-7,1.0e-7,100)

%3-2-1
 %P1= 1.192309125394795   1.221601050222906
 %P2=  -0.286032169310697  -0.118185598207315
 
 %3-2-2
 %P=-0.269503233982088  -1.315483472346922
 %P= 0.669951470894617   1.633998614104253